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ABSTRACT 

We propose a simple model for the origin of fast and slow rotator early-type galaxies 
(ETG) within the hierarchical ACDM scenario, that is based on the assumption that 
the mass fraction of stellar discs in ETGs is a proxy for the specific angular momentum 
expressed via \r- Within our model we reproduce the fraction of fast and slow rotators 
as a function of magnitude in the ATLAS '^^ survey, assuming that fast rotating ETGs 
have at least 10% of their total stellar mass in a disc component. In agreement with 
ATLAS'^^ observations we find that slow rotators are predominantly galaxies with 
M* > 10^°-^ M0 contributing ~ 20% to the overall ETG population. We show in 
detail that the growth histories of fast and slow rotators are different, supporting 
the classification of ETGs into these two categories. Slow rotators accrete between 
~ 50% — 90% of their stellar mass from satellites and their most massive progenitors 
have on average up to 3 major mergers during their evolution. Fast rotators in contrast, 
accrete less than 50% and have on average less than one major merger in their past. 

We find that the underlying physical reason for the different growth histories is the 
slowing down and ultimately complete shut-down of gas cooling in massive galaxies. 
Once cooling and associated star formation in disc stops, galaxies grow via infall from 
satellites. Frequent minor mergers thereby, destroy existing stellar discs via violent 
relaxation and also tend to lower the specific angular momentum of the main stellar 
body, lowering Xr into the slow rotator regime. 

On average the last gas-rich major merger interaction in slow rotators happens 
at z > 1.5, followed by a series of minor mergers. These results support the idea that 
kinematically decoupled cores (KDC) form during gas-rich major mergers at high-z 
followed by minor mergers, which build-up the outer layers of the remnant, and make 
remnants that are initially too fiat compared to observations become rounder. Fast 
rotators are less likely to form such KDCs due to the fact that they have on average 
less than one major merger in their past. 

Fast rotators in our model have different formation paths. The majority, 78% has 
bulge-to-total stellar mass ratios B/T > 0.5 and managed to grow stellar discs due to 
continued gas cooling or bulges due to frequent minor mergers. The remaining 22% 
live in high density environments and consist of low B/T galaxies with gas fractions 
below 15%, that have exhausted their cold gas reservoir and have no hot halo from 
which gas can cool. These fast rotators most likely resemble the flattened disc-like fast 
rotators in the ATLAS'^^ survey. 

Our results predict that ETGs can change their state from fast to slow rotator 
and vice versa, while the former is taking place predominantly at low z (z < 2), the 
latter is occurring during cosmic epochs when cooling times are short and galaxies 
gas-rich. We predict that the ratio of the number density of slow to fast rotators is a 
strong function of rcdshift, with massive (> 10^*^ M©) fast rotators being more than 
one order of magnitude more frequent at z '^ 2. 

Key words: galaxies: elliptical and lenticular, cD - galaxies: formation - galaxies: 
evolution - galaxies: structure 



1 INTRODUCTION 



Much attention has been paid to the modelling of 
the formation of early- type galaxies (ETG) o v er the 
last th ree decades (e.g. pToomre fc Toomrel Il972l: IWhite 



19781 : iNegroponte fc Whitd 1 19831: iBarnes fc Hernguist 



1992: iMihos fc Hernauistl Il996l: JNaab fc BurkertI 1200^ 



Bournaud et all 12004 : ICox et al.l l2006l : iHoffman et all 
2010). Once believed to be coeval featureless stel- 
lar systems in virial equilib r ium that formed from a 
'monolithic' collapse l|Larsonl 11974 ). an ever- increasing 
wealth of observational data has revealed a multitude 
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of distinct physical properties: in the ir rotational sup- 
port ( Binney 197a iDavies et al.l 19831). i s ophot al shape 



metric profiles ( Ferrarese et al 



iKormendv et al 



iBender et all 1 19891: iKormendv fc Bender! Il996l), photo 



v &g aenaea iiijijqi. pnoto- 
1.1 Il994l : iLauer et al.l Il995l : 



20091). and in t heir stellar angu l ar mo - 



mentum (JEmsellem et all l2007l : ICappellari et all |2007| ). 
Driven by such observations theoretical models have been 
put to the test, and in particular the ACDM paradigm 
for structure formation itself within which structure grows 
hierarc hically via a seque nce of mergers and accretion 
events (|White fc Reeslll978l ). 

Early nu merica l N-body simulations 

iToomre fc Toomrel 1 19721 ) identified galaxy mergers as 
a viable mechanism to transform dynamical cold stellar 
systems into dynamical hot ones resembling early-type 
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galaxies with a de Vaucou l eurs-like light pr o file (s ee also 
iBarnes fc HernquistI Il992l : iNaab fc Truiilld |2006| ). The 
main driver for this^ morphological transformation being 
violent relaxation ( |LYndcn- Bcll 1967) during which stars 
get redistributed. The merger scenario subsequently has 
been tested against two main categories of observational 
data: global propert ies in large homogenous surveys (e.g. 
iBernardi et al.ll2003l ) and detailed high-resolution observa- 
tions of an inhomog eneous set of individual galaxies (e.g. 
Ide Zeeuw et al.ll2002l v 

Galaxy formation models, in particular semi-analytic 
models, employing the merger scenario, reproduce success- 
fully the mass function, colour distribution, star formation 
history and metallicity of early-type galax ies as observed 
e.g. i n the Sloan Digital Sky Survey (SDSS) (|De Lucia et al.l 
[200g, e.g.). However, they still have problems in e.g. re- 
coveri ng the trend in a/Fe with stellar velocity disper- 
sion (JThomas et al.1 l2005l : iNagashima et al.1 120051 ') or the 
j oint distribution in size s and velocity dispers ion at z = 
l|van der Wei et al.1 [2OO9I : IShankar et al.1 I2OI0I ). The latter 
has also been reflected in mismatches to the fundamental 
plane of early type galaxies d Almeida et al. 2 007 ). While the 
deviations regarding the sizes and velocity dispersions likely 
stem from too efhcient in-situ star f ormation in progen i- 
tors of present-day early- type galaxies IjShankar et al.ll2010l ) , 
the failure to produce enough a-elements is most likely con- 
nected to a combinat ion of too much residu al star formation 
dThoinas et al.ll2005h and satellite infall iJKhochfar fc Silkl 
l2006al ). 

The inner parts of early-types show further detailed 
characteristics in t he form of steep increasing power-law 
or core like profiles IJFerrarese et al.lll994l : iLauer et al.ll 19951 : 
iKormendv et al.ll2009l ). The former has been associated with 
cold gas in discs that loses its angular momentum during the 
merging process (e.g. jBarncs fc Hornquist 1991) and set- 
tles to the centre of the remnant where it triggers a star 
burst, thereby steepening the light profile (jHopkins et al.l 
[2003). Core galaxies on the other hand are most likely 
the product of a non-dissipative merger event between two 
galaxies hosting super-massive black holes (SMBH) at their 
centr e (iMilosavlievic fc MerrittI I2OOII : iKhochfar fc BurkertI 
I2OO5I : iHopkins et al.l l2009d ). During the inward spiralling 
the SMBH pair kick out stars via three-body interaction 
and p roduce s a density core in the centre of the remnant 
l|Milosavlievic fc Merrittll200ll '). 

Classifying elliptical galaxies by the deviation of their 
isophotal shape from pure ellipses has been suggested as 
a method to separate the population of elliptical galaxies 
into two classe s with box-like (b oxy) and disk-like devi- 
ations (disky) IjBender et al.lll989l V Numerical simulations 
show that equal mass mergers in general produce remnants 
with boxy isophotes due to stronger gravitation al inter- 
action and associated violen t relaxation (Bcndo fc Barnea 
I2OOOI : iNaab fc BurkertI bOOJ ) efficiently destroying progen- 
itor discs. In contrast unequal mass major mergers with 
2 ^ M1/M2 ^ 4 result on average in disky remnants 
for the majority of projectio ns under which they are 
viewed (jNaab fc BurkertI 120031 ). Observationally, the aver- 
age isophotal shape of elliptical g alaxies turns toward b oxy 
going up in mass (see howerver lEmsellem et al.l I2OIII ) in- 
dicating a transition in the formation process of massive 
ellipticals that cannot be attributed to the mass ratio dur- 



ing the last major merger alone, which in general does not 
chang e significantly enough in frequency as a function of 
mass (|Khochfar fc B urkert 20051). In addition unequal mass 
mergers of early-type galaxies are able to transform disky 
ellipticals into box y ones in the majority of simulated cases 
(JNaab et al.ll2006r ). Within a cosmological context the rela- 
tive fraction of dry, early-ty pe major mergers increases with 
stellar mass of the remnant (Khochfar fc Burkert 2003) pro- 
viding the additio nal channel needed to recov er the observed 
trend with mass JKhochfar fc BurkertI |2005| ). In a detailed 
study of the orbital content of stars in merger remnants 
Jes seit et al.l (|2005l ) showed that the individual mixture of 
orbits is the deciding factor on the isophotal shape of disky 
and boxy elliptical galaxies. 

Making full use of the dynamical information avail- 
able in integral-field observations the SAURON sample 
Jde Zeeuw et al.l 120021 ) measured the specific baryonic an- 
gular momentum content via the proxy Xa within one Re 
of a representative sample of early-type galaxies and sepa- 
rate them into fast and slow rotators according to their A_r 
value, in a r obust way that i s nea r ly i nsensitive to projec- 
tion effects dEmsellem et al.l [2OO7I : ICappcUari ct al. 200'(j). 
Observationally, fast rotators present regular rotation pat- 
terns aligned with the photometry, while slow rotators have 
low angular momentum content and show misalignments 
between the photometry and the velocity axes and ofte n 
exhibit kinematically distinct cores (JKrainovic et al.ll2008l ). 
The division into fast and slow rotators is fur ther supported 
by recent simulations of (|Bois et al.l l2011al . hereafter Pa- 
per VI) that show that fast (resp. slow) rotators have a 
high (resp. low) angular momentum content, their photo- 
metric and kinemetric position angle are aligned (resp. mis- 
aligned), they present (resp. do not present) regular veloc- 
ity patterns. W hile the SAURON surv ey was biased towards 
massive ETGs (Ide Zeeuw et al.ll2o"o3 ) the ATLAS^°survey 
ICappellari et al.ll2011al . hereafter Paper I) for the first time 
allows to give a complete view of the local ETG popula- 
tion within a radius of 42 Mpc using the full information of 
integral- field observations. The combination of completeness 
with detailed structural info rmation and comp lementary ob- 
servations on e.g. the HI ijSerra et al.ll201ll ) or molecular 
I Young et al.ll201ll ) gas masses, make the ATLAS ^"^ survey 
an ideal data set to investigate the formation of ETGs. 

In this paper we attempt to model fast and slow rotators 
within the standard ACDM paradigm and analyse specific 
differences in their formation. The structure of the paper is 
as follows: we start by summarising the basic model ingre- 
dients used in our semi-analytic modelling (SAM) approach 
in section O followed by showing the performance of the 
model in section [3| In section |4] we introduce our model for 
fast and slow rotators and investigate the robustness of our 
assumptions before we analyse in detail differences in their 
formation in section [5] and predict their redshift evolution in 
section[Sl Section[7]is devoted to the question of cosmic vari- 
ance in the ATLAS'^^sample followed by sections [S] fc [5] in 
which we summarise our results and conclude respectively. 



2 THE MODELLING APPROACH 

We apply the semi-analytic modelling (SAM) approac h to 
galaxy formation pioneered in a series of papers (e.g. ISilkl 
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ll977l : IWhite fc ReedflQTJ : IWhite fc Frenklll99lh . The main 
strategy behind SAMs is to first calculate the collapse and 
merging history of individual dark matter halos, and sec- 
ondly to approximate the more complex physics of the 
baryons inside these dark matter halos by analytic models 
(see e.g-lKauffmann et alJI 19991: Somerville fc Primackl 199E : 



Cole et all I2OO0I: iHatton et all 20031: ICattaneo et all 



I 



Croton et alj '2006: Bower et al." "2006"; "Kang _et all l200e ; 
De Lucia fc Blaizot 2007 : Khochfar fc Silk 2009bl ). Each of 



2ooe 



the dark matter halos will consist of three main components 
which are distributed among individual galaxies inside them: 
a stellar, a cold, and a hot gas component, where the latter 
is only attributed to central galaxies, which are the most 
massive galaxies inside individual halos and are typically 
obse rved to reside in extended X -ray emitting hot haloes of 
gas IjMulchaev fc Jeltemall201Gl ). In the following sections, 
we will describe briefly the recipes used to calculate these 
different com ponents which are mainly based on recipes pre- 
sented in e.g. iKauffmann et alJ l| 19991 ) (hereafter, K99) and 
IStaringel et al.l (200l|) (hereafter, SOI), and we refer read- 
ers for more details on the basic model implementations to 
their work and references therein. In the remainder of this 
paper we call SAM the model im plementation presented in 
iKhochfar fc Burkerd l|2003l, 120051 ) which is summarised in 
the following parts of this section. "With our SAM we model 
a cosmic volume of Va3d = 10^ Mpc^ , which is comparable 
to the volume of the ATLAS^°survey (Paper I). Throughout 
this paper, we use the following set of cosm ological param- 
eters consistent with 7-year WMAP data l|Komatsu et al.l 
I2OIOI ): fio = 0.26, ^A = 0.74, ^5/^0 = 0.17, Us = 0.97, 
o-g = 0.8 and h = 0.71. 



2.1 Dark Matter 

The mass function of dark matter haloes within our vol- 
ume is calcu l ated using the analytic fitting formula of 
Ijenkins et al.l l|200ll ). In addition, to investigate the effects 
of cosmic variance (see Section [7| we have run a cosmo- 
logical N-body simulation of a (100 Mpc//i)'' volume with 
GADGET-2 (Springel 2005 ). The particle resolution of the 
simulation is 10^ Mq/Zi, and we generate the z = Q halo 
mass function using the friends-of-friends algorithm. 

The corresponding halo merging histories are calculated 
using a Monte-Carlo appr oach based on the Ex tended-Press- 



Schechter formalism (e.g. Lacev fc Cold 19931) fol lowing the 
method presented in ISomerville fc KolattI ( 19991 ) . This ap- 



proach has been shown to produce merging histories and 
progenitor distributions in reasonable agreement with re- 
sults from N-body simulations of col d dark matter structur e 
formation in a cosmological context IjSomerville et al.ll2000^ . 
The merging history of dark matter halos is reconstructed 
by breaking each halo up into progenitors above a limiting 
minimum progenitor mass Mmin- This mass cut needs to be 
chosen carefully as it ensures that the right galaxy popula- 
tion and merging histories are produced within the model. 
Progenitor halos with masses below Mmin are declared as 
accretion events and their histories are not followed further 
back in time. Progenitors labelled as accretion events should 
ideally not host any significant galaxies in them and be com- 
posed mainly of primordial gas. To achieve a good compro- 
mise between accuracy and computational time, we estimate 
Mmin by running several simulations with different resolu- 



tions and chose the resolution for which results in the galaxy 
mass range of interest are independent of the specific choice 
of Mmin- Changing the mass resolution mainly affects our 
results at low galaxy mass scales, leaving massive galaxies, 
the focus of this work, literally unaffected. Throughout this 
paper we will use Mmin = 10^° M© which produces numer- 
ically stable results for galaxies with stellar masses greater 
~ 6 X 10^ M0, which corresponds to the mass limit in the 
ATLAS'^^survey. For simplicity we assume in the following 
that the dar k matter profiles are th at of truncated isother- 
mal spheres (jKauffmann et al.lll999f ). 



2.2 Baryonic Physics 

Once the merging history of the dark matter component has 
been calculated, it is possible to follow the evolution of the 
baryonic content in these halos forward in time. We assume 
each halo consists of three components: hot gas (Mhot), cold 
gas (Mgas) and stars (Af»), where the latter two compo- 
nents can be distributed among individual galaxies, labelled 
central and satellite, inside a single dark matter halo. The 
hot gas component however, is only associated with central 
galaxies. The stellar component of each galaxy is addition- 
ally divided into bulge (Mtui) and disc {Mdiac, to allow mor- 
phological classification of model galaxies. In the following, 
we describe how the evolution of each component is calcu- 
lated. 



2.2.1 Gas Cooling & Reiomsation 

Each branch of the merger tree starts at a progenitor mass 
of Mmin and merges at a redshift 2^0 with the main 
progenitor branch. Initially, each halo is occupied by hot 
primordial gas which is captured in the potential well of 
the halo and shock heated to its virial temperature Tvir = 
35.9 rV^/ (km s~'^)1 K, where Vc is the circular velocity of 
the halo (iwhite fc Frenklll99ll . K99). Subsequently this hot 
gas component is allowed to radiatively cool and settles 
down into a rotationally supported gas disc at the centre 
of the ha l o, which we identify as the central galax y (e.g 
[Siii3ll977l : IWhite fc Reeslll97a : IWhite fc Frenk|[l99ll ). The 
rate at which hot gas cools down is estimated by calcu- 
lating the cooling radius inside the halo usin g the cool- 
ing functions provided bv [Sutherland fc Dopital (119931 ) and 
the prescription in SOI assuming the gas density profile fol- 
lows that of the hosting dark matter at all times. At early 
times and in low mass haloes the cooling times of the gas 
are much shorter than the halo dynamical time suggesting 
that no stable shock will exist at the virial radius of the 
halo, and th at the gas d ocs not shock-heat to virial tem- 
perature ( Birnboim fc Dc kcl 2003; Kcrcs ct al. 2005;). We 
capture this situation by calculating the radiative cooling 
time of the halo gas under the assumption of it being at 
virial temperature and comparing it to the halo dynami- 
cal time. Whichever is longer will be used to calculate the 
cooling rate/accretion rate of gas onto the central galaxy. 
This procedure has been shown to provide an accurate dis- 
tinction between haloes with shock heat ed gas and without 
around a halo mass scale o f ~ 10^^ M© (jCroton et al.ll2006l : 
IKhochfar fc Qstrikeill2008l ). Haloes more massive than this 
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Figure 1. The distribution of SAM accretion rates of cold gas 
onto central galaxies living in halos with MoAi ~ 10^^ Mq at 
z = 2.5. The distribution is in good agreem ent with the results 
from fully numerical simulations presented in lPekel et al.l ll2009|). 
Please note that we do not include major mergers in this plot due 
to the ambiguity of the time scale over which the gaseous mate- 
rial of the satellite galaxy contributes to the central galaxy. The 
missing m ergers expla i ns the slight mismatch at high accretion 
rates with lDekel et al.l l l200g|) . 



show at z ^ 2 signs of accretion of cold gas al ong cosmic fil- 
aments (jOcvirk et al.ll2008l : IPekel et al■ll2009^ ■ which is not 
captured within the a bove model ling approach. We here in- 
clude the approach of lKhochfar fc Silk (2009b) that models 
this behaviour and is based on f i ts to t he results from numer- 
ical simulations bv lKeres et al.l (|2009|). In Fig.[l]we show the 
distribution of accretion rates onto central galaxies in a halo 
of Mdm ~ 10^^ M© at z = 2.5. The distribution is in good 
agreement w ith the resul t s from fully numerical simulations 
presented in lDekel et al.l (|2009l '). 

In the case of a merger between halos we assume that 
all of the hot gas present in the progenitors is shock heated 
to the virial temperature of the remnant halo, and that gas 
can only cool down onto the new central galaxy which is 
the central galaxy of the most massive progenitor halo. The 
central galaxy of the less massive halo will become a satel- 
lite galaxy orbiting inside the remnant halo. In this way, a 
halo can host multiple satellite galaxies, depending on the 
merging history of the halo, but will always only host one 
central galaxy onto which gas can cool. The cold gas content 
in satellite galaxies is given by the amount present when 
they first became satellite galaxies and does not increase, 
instead it decreases due to ongoing st ar formation and su- 
pernova feedback (see however, e.g . iKhochfar fc Ostrikeii 
120081 : iKang fc van den BoschI 120081 . for models on gradual 
gas stripping). 

In the simplified picture adopted above, the amount of 
gas available to cool down is only limited by the universal 
baryon fraction. However, in the presence of a photoionising 
backgroun d the fraction of baryons captured in small halos 
is reduced l|GnedirJl2000l ) and we use the recipe of lSomervillel 



(|2002r ) , which is based on a fi tting formulae derived from hy- 
drodynamical simulations bv lGnedirJ ((20001), to estimate the 
amount of baryons in each halo. For the epoch of reionisa- 



tion, we assume Zr 



10, based on recent results from 



WMAP JKomatsu et al.ll2010h . 



2.2.2 Cooling shut-off mass scale 

Cooling rates for massive galaxies must drop off significantly 
in order to e.g. reproduce the obs erved colour distrib ution 
and mass function of galaxies (e.g. ICroton et al.ll2006l ) . Pos- 
sible candidates for physical processes responsible in shut- 
ting down cooling in massive dark matter haloes are e.g . 
AGN- feedback (e. g lGranato et al.ll2004l : ICroton et a"i]|2006l : 
I Bower et al. 120061) or gravitational heat i ng by infalling sub- 
structures ( Kho chfar fc Ostriker 2003; iDekel fc BirnboimI 
120081 : [Johansson et al. 2009 ). Effectively these processes re- 
sult in cooling being suppressed in haloes above a critical 
mass scale, and have been implemented using a cooling shut- 
off in haloes above a critica l circular velocity (K99) or mass 
scale (jCattaneo et al.ll200q ). By applying an effective shut- 
off mass scale of MnM,crit ~ 10^^ M0 at z ^ 2 one is able 
to repro duce the colour distrib ution and mass function of 
galaxies IjCattaneo et al.ll200q ). Additionally, such models 
are able to r eproduce the frequency of dry, gas-poor, ma- 
jor mergers f Kh ochfar fc Silkll2009al ). In the following we 
will modify our fiducial cooling model as laid out in chap- 
ter [2]2Tl] by adopting a quenching of cooling in dark matter 
haloes above a critical mass scale of MDM,crit ^ 2.5 x 10^^ 
Mb at z ^ 2 which is similar to the value proposed in 
ICattaneo erall (|2006h {MDM,crit ^ 2 x 10^^ Mq), and al- 
lows to have gas in massive late- type galaxies (see Fig. U)). 
Note that we allow the gas that is already in the disc of a 
galaxy to continue forming stars until it is used up, even 
after the host halo is more massive than Mdm cru. 



2.2.3 Star Formation in Discs and Supernova Feedback 

Once cold gas has settled down in a disc, we allow for frag- 
mentation and subsequent star formation accordi ng to a pa- 
rame terised global Schmidt-Kennicutt type law (JKennicuttl 
[1998) of the form M* = aMgas/tdyn,gai, where a is a free 
parameter describing the efficiency of the conversion of cold 
disc gas Mgas into stars, and tdyn,gai is the dynamical time of 
the galaxy and is approximated to be 1/10 ti mes the dynam- 
ical time of the hosting dark matter halo ||Mo et 



i d ynam- 
aL[ |l998l . 



K99) . For a satellite galaxy we calculate the dynamical time 
using the halo properties when it last was a central galaxy. 
We assume here that newly formed stars M« all contribute 
to the stellar disc Mdisc of the galaxy. 

Feedback from supernovae plays an important role in 
regulating star formation in small mass halos a nd in prevent- 
ing too massive satellite gal axies from forming (JDekel fc Silkl 
1 19861 : iKhochfar et~aLll20()i1 ). We implement feedback based 
on the prescription presented in K99 with 

AM....., = 1.^Z£^AM,. (1) 

O Vc 

Here we introduce a second free parameter e which repre- 
sents our lack of knowledge on the efficiency with which the 
energy from supernovae reheats the cold gas. The expected 
number of supernovae per solar mass of stars formed for a 
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typical IMF (e.g. Scalo) is t^sn = 5 x 10"^, and Esn = 10^^ 
erg is the energy output from each supernova. We take Vc 
as the circular velocity of the halo in which the galaxy was 
last present as a central galaxy or is right now. 

2.2.4 Implementation 

In between mergers the evolution of galaxies is governed by 
the equations presented in the previous sections. We solve 
this coupled system of differential equations by integration 
in between time steps of constant Az = 0.02, during which 
the halo properties are kept fixed. We thus solve following 
set of equations: 



M. 



M„ 



^dyn,gal 



Mreheat 


= 


3' K^ M, 


Mgas 


= 


Mcool — Mreheat " 


Mhot 


= 


— Meool + Mreheat 


Mdisc 


= 


M,. 



(2) 



M, 



Here Mcooi is the cooling rate applying the prescriptions 
of the previous sections. The above set of equations is strictly 
only applicable to central galaxies, which show cooling of hot 
gas. For satellite galaxies we use a modified set of equations 
with Mhot ~ Mhot ~ Meooi ~ 0. All gas that is reheated 
in satellites is added to the hot gas reservoir of the central 
galaxy. 

2.2.5 Galaxy Mergers and Bulge Formation 

We allow for mergers between galaxies residing in a single 
halo. As mentioned earlier, each halo is occupied by one cen- 
tral galaxy and a number of satellite galaxies depending on 
the past merging history of the halo. Whenever two halos 
merge, the galaxies inside them will merge on a time-scale 
which we calculate by estimating the time it would take 
the satellite to reach the centre of the halo under the effect 
of dynamical friction. Satellites are assumed to merge only 
with central galaxies and we set up their orbits in the halo 
according to the prescription of K99, modified to use the 
Coulomb log arithm approximat ion of SOI and the correc- 
tion factor of I Jiang et al.l ([2008'). Expressed in units of the 
Hubble time the dynamical friction is only dependent on the 
mass ratio of the merging partners. The merger rates in our 
SAM are in good agreement with observations (jJogee et al.l 
|2009|'1 and agr ee within a factor of few with other models 
l|Hopkins et al.ll2010al lbtl. 

We do not distinguish between the effects of major 
R = Mgai,i/Mgai,2 < 3.5 {Mgai,i > Mgaia) and minor 
mergers 10 > -R > 3.5 in our SAM, but rather apply a 
model that shows a smooth transition between the effects 
of such mergers. In what follows the index 1 will always 
refer to the primary galaxy and the index 2 to the merg- 
ing satellite galaxy. During major mergers the stellar discs 
of the progenitors will be almost completely destroyed and 
added to the bulge component of the remnant. In practice 
we assume that the amount of disc mass scattere d into the 
bulge is mm[Mdisc,i,M,^2] (|Hopkins et al.ll2009al 'l. In addi- 
tion we assume that a star burst will be triggered due to 
the interaction of the galaxies. The fraction of cold disc gas 



transformed into stars during the star burst is modelled us- 
ing a function fturst that includes dependencies on the cold 
gas fraction fgas = Mgas /{Mgas + Mt), the disc fraction 
fdisc = Mdisc I M* and mass ratio R of the merging galaxies 
iHopkins et al.ll2009al ^ . We here use the following functional 
form that shows broad agreem ent with results from d etailed 
numerical merger simulations (jHopkins et al.ll2009al ). 



J bur 



1 — (1 + rsb/rdisc) exp{—rsb/rdisc). 



(3) 



We calculat e the d isc scale length rdisc using the disc model 
of lMo et al.l (|l998l ) and the radius rsb, within which gas loses 
its angular momentum and contributes to the central star 
burst via: 



rsb 



"s6(l — fgas)fdiscFoG{R), 



(4) 



with G{R) = 2R/{1 + R). We set the numerical con- 
stant asb ~ 1 and the orbit dependent constant Fg — 1.2 
iHopkins et al.l l2009al ). The latter represents an average 
value for a range of different merging orbits. 

The cold gas that is not converted into stars during the 
merger is added to the gaseous disc of the remnant. This 
results in following composition for the remnants of major 
and minor mergers: 



M, 



bul.rem — 



■^* -^ disc, rem- 
■^''J- g a s ,r em 



Mbui.i + Mbui,2 + mm[Mdisc,i,M^,2] 

+ Mdisc,2 + fburst,lMgas,l 

+ fburst,2Mgas,2 (5) 

Mdisc.i - mm[Mdisc,i,M,,2] (6) 

(1 ~ .fburst,l){Mgas,l) 

+ (1 - .fburst,2)Mgas,2). (7) 



Immediately after major mergers our model galaxies 
have a large fraction of their stars in the bulge component 
and the stellar disc can start re-growing by either conversion 
of cold gas that was not used up in the merger or by freshly 
accreted gas. We here assume mergers as the only efficient 
means of bulge formation and neglect other possible secula r 
effects such as e.g. disc instabilities (see e.g. ICombesll2009l '). 



3 MODEL NORMALIZATION 

The various prescriptions for physical process introduced in 
the last section include a number of free parameters which 
represent our lack of knowledge on parts of the underlying 
physical processes. Generally SAMs use a set of observations 
to fix these parameters. Reproduction of such observations 
thus cannot been seen as a model prediction but rather a 
'tuning' of it. In this section we present the core observations 
our free parameters are tuned to and their individual values. 



3.1 Mass function 

The mass function of galaxies is a core observable that is well 
measured and robust. Much progress in galaxy modelling 
has been achieved over the last years by focusing on repro- 
ducing specific features in it such as e.g. the s t eep exponen- 
tial decline at high masses (IBenson et al.ll2003l : ICroton et al.l 
l2006l : iKhochfar fc Ostrikerl bOOSJ ) or the s teepness of the 
faint- end slope as a function of redshift (jKhochfar et al.l 
l2007f ). The ATLAS^°parent sample (Paper I) shows very 
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Figure 2. The mass function of galaxies at 2: = in our best-fit 
model (solid line). Filled circles are measurements by I Bell et al.l 
l|2003l V Red and blue lines show the early- type and late- type 
galaxy mass function, respectively. 



good agreement with the mass function of lBell et al.l ((2003), 
to which we compare our best-fit z = model mass function 
(Fig. [2]). We find very good agreement both in shape and 
overall normalisation over a wide range in masses, in partic- 
ular at the high mass end, where our model population is 
dominated by early-type galaxies with bulge-to-total mass 
ratios B/T > 0.7 (see Fig. [3]). 







JS^sf*****'*^ 



9.0 9.5 10.0 10.5 11.0 11.5 12.0 
log(MyMQ) 

Figure 3. Bulge-to-total mass ratios B/T of galaxies as a func- 
tion of mass. The solid blue line shows the median value at a 
given galaxy mass. Our results show a clear trend towards bulge- 
dominated galaxies at larger masses. The overall scatter of the 
distribution is largest around 10"'^"'^ M0, which also marks the 
point at which the galaxy populations is largely dominated by 
ea rly- type galaxies. T he solid triangle shows the Milky- Way value 
of lSmith et al.l 1 120071 ) . The horizontal solid and dashed line mark 
the transition to early-type galaxies {f early = 0.5) and slow ro- 
tators {/fast = 0.1), respectively. Large red circles show the final 
extended early-type galaxy sample in the SAM, including the 



additional criterion f^e 



Mgas/M, = 0.15, on the minimum 



gas fraction in late- type galaxies. The latter has been included to 
mimic the ATLAS^^ selection of early-type galaxies, which allows 
for the inclusion of gas-poor low B/T galaxies. 



3.2 BuIge-to-Total Mass Ratios 

Continued merging is transforming morphologies of galaxies 
by scattering disc stars into bulges. If modelled correctly, 
the observed trend of increasing bulge fractions with galaxy 
mass should be recovered by the model. Since we focus here 
on early-type galaxies this is an important test. In Fig. [3] 
we show the B/T and stellar masses of our model galaxies, 
as well as the median B/T at a given mass. We are able 
to recover the general observed trend toward bulge dom- 
inated systems with increasing mass. The median B/T is 
below 0.2 for galaxies with A'l, ^ 5 x 10^" M© and makes 
a strong transition to B/T > 0.6 at M, > lO" Mq. It ap- 
pears when the cooling of gas starts slowing down, in haloes 
below MoM.crit- As a consequence discs grow slower, while 
merging continues to transform discs into spheroids, shifting 
the median B/T to larger values. 



3.3 Gas Fraction 

Another important quantity related to the efficiency of disc 
re-growth and survival during mergers is the cold gas frac- 
tion, fgas = M gas /{M gas + M-j,), as a function of the stellar 



mass. We show the median value from our model for late- 
type galaxies with B/T < 0.5 in Fig. [H and compare it to 
the observed range in disc galaxi es (JBell fc de Jond I2OO1I : 



c ga , 

iKannapparil I2OO4I : iMcGaughl I2OO5I ) . We find fair agreement 



for model galaxies with M* > 10 M© and the observations. 
The decreasing gas fraction in our model is mainly due to 
the increasing cooling time as a function of host halo mass, 
and the sharp cooling shut-off in haloes above MoM.crit- 



3.4 Early-Type classification 

In general SAMs separate between early-type and late- 
type galaxies using the bulge-to-total ratio B/T either 
in nia,ss or _B— band ma g nitude (e.g. iKauffmann et al.l 
I1999I : ISpringel et al] I2001I : iBower et al] bood ). This ap- 
proach is based on the good correlation between the 
Hubble T-type and the i3— b and B/T ratio of galaxies 
(jSimien fc de Vaucouleurslll986l ) . The early- type galaxies se- 
lected within the ATLAS^'' survey are in very good agree- 
ment with morphologies based on the T-type classification 
scheme (Paper I). We will use here a fiducial value for the 
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bulge-to-total mass ratio of B /T — 0.5 to separate early- 
type from late-type galaxies. 

The ATLAS^'^survey however, contains a significant 
numbers of flat disc-like galaxies, that have morphologies 
similar to late-type spirals, but do not show evidence for 
gas, dust or spiral arms, and are for this reason classified as 
early-types (|Cappellari et al.ll2011b|) according to the stan- 
dard morphological classification ( Sandage 1 1961 ). Such red 
spirals are classified a s fast rotators ( Krainovic et al.ll201ll : 
lEmsellem et al.ll201ll . hereafter Paper II and Paper III, re- 
spectively) and are not included, if the early-type selection 
in the SAM is done purely based on a B/T cut. The com- 
bination of a strong decline in the mass function of galaxies 
and the transition from fast to slow rotator dominance, sug- 
gest that the impact on the ratio of fast to slow rotators by 
missing such galaxies, could be important at masses ^ 10^^ 
M© (Paper III). We address this issue by including an ad- 
ditional criterion based on the gas fraction in galaxies with 
low B/T. We first compare the distribution of the cold HI 
gas fraction fred = Mgas/Mt. in model ce ntral galaxies with 
that obtained from the ATLAS^°sample (|Serra et aLlbOllh 
assuming a constant mo lecular-to-atomic ratio (see however, 
lObreschkow et al.ll2009l . for more detailed modelling of the 
molecular fraction in cold gas). The distributions show very 
similar trends with a peak at high gas fractions for late-type 
galaxies and ETGs spanning a wide range from low to high 
gas fractions with a strong tendency for low gas fractions 
(Serra et al.ll201ll '). Modelled late-type galaxies {B/T < 0.5) 
show in addition a small peak at gas fractions below ~ 15% 
which is not seen in the observations Such galaxies are likely 
candidates for fast rotators with low B/T and we therefore 
adopt an additional criterion for early-type galaxies with 
B/T < 0.5 and fred < 0.15. A significant fraction of galax- 
ies with masses ~ lO^^M© in our model are centrals liv- 
ing in haloes which have been crossing MDM,crit some time 
ago, and since then did not encounter any significant merg- 
ers thus using up their gas reservoir and building a stel- 
lar disc with little interference. The hosting halo of these 
galaxies are not very much larger than Mdm^ctu, thus ex- 
plaining why they do not reach masses larger than ~ 10^^ 
M0. In general we find that the majority of low B/T fast 
rotators below ~ lO^^M© are satellite galaxies. The frac- 
tion of early-type galaxies is known to inc rease with the 
density of the environment (Dressier 1980). A consistent 
trend has been found for the galaxies in the ATLAS'^^survey 
ICappellari et al.ll2011bl . hereafter Paper VII). To include to 
first order gas-poor low B/T satellite galaxies in high density 
environments, we add satellite galaxies in group-like dark 
matter haloes with Mdm ^ 4 x 10^"^ M© and stellar masses 
larger than 10^" M©. Again we here use the early- type def- 
inition with B/T < 0.5 and fred < 0.15 for these satellite 
galaxies. We find that the trend in the ratio of fast to slow 
rotators is not significantly affected by adding the satellite 
population, as most of them have masses in a range where 
fast rotators clearly dominate the population of early-type 
galaxies. 

In Fig. Owe show the final selection of early-type galax- 
ies from the overall model sample. This selection yields a to- 
tal number density of 1.6 x 10~^Mpc~^ early-type galaxies 
with M, J5 6 X 10^ Mq, which is in good agreement with 
the observed value for the ATLAS^°sample of 2.2 x 10"^ 
Mpc~^ (Paper I). In Fig. [2]we show the mass function sepa- 



star formation efficiency 
supernovae feedback efficiency 
minimum gas fraction in late-types 
early-type threshold 
fast rotator threshold 



Jred 
Jearly 
J fast 



0.02 

0.4 

0.15 

0.5 

0.1 



5 [33] 
3 [Ml 



Table 1. Most important free model parameters and their 
adopted values. 
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Figure 4. The median fraction of gas fgas = Mgas/{Mgas +Mt) 
as a function of galaxy mass in the SAM (filled circles) range cov- 
ered by observations including their error bars (hatched region). 
Solid squares show the compilation of various ob servational data 
JBell fc de Jonell200ll : lKannappa32004l : lMcGaugh 2005 ). and the 
solid triangle is the value for the Milky- Way. The declining gas 
fraction with stellar mass is mainly due to the increasing cooling 
time scales for the gas. Please note, that due to the sharp cooling 
shut-off in halos above MoM,crit the gas fraction declines very 
sharply to zero. 



rated into early-type and late-type galaxies for our adopted 
selection. The high-mass end is as expected dominated by 
ETGs, and we find that the number density is nearly con- 
stant for M, < 10^^ Mq. Compared to the observed mass 
function presented in Paper I our model produces too few 
ETGs at low masses, most of which should be fast rotators 
and too few very massive M« > 10^^ M© late-type galaxies. 
Approximately 80% of the model early-type galaxies 
have B/T > 0.5 while the remaining 20% are low B/T 
galaxies devoid of gas. 



4 FAST AND SLOW ROTATING ELLIPTICAL 
GALAXIES 

The results of the ATLAS^^survey suggest that ETGs can 
be classified into two classes based on the parameter \r = 
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{R\V\)/{Ry/V"^ + a'^) which has been shown to be a good 
proxy for the specific baryonic angular m omentum in ETG 
IJEmsellem et al.ll2007l : ljesseit et al.ll2009l . Paper III). Using 
a suite of high resolution simulations Paper VI shows that 
during major mergers (R < 2) of late-type galaxies the ge- 
ometry of the encounter is the dominant factor in deciding 
whether the remnant is a slow or fast rotator. In general, 
major mergers that are retrograde (prograde) with respect 
to the main progenitor result in slow (fast) rotators. In mi- 
nor mergers (R Js 3) remnants are almost exclusively fast 
rotators. However, the situation is not as simple as that. 
Simulated slow rotators tend in general to be too flat com- 
pared to observed ones and have properties similar to 2(t 
galaxies (Paper II), suggesting that a single major merger 
might not be the way to form them. Furthermore, merg- 
ing of remnants show that the orbital angular momentum is 
able to spin up slow rotators, even in major mergers. Thus 
to predict the proper outcome of a binary merger, one needs 
to know the initial state of the progenitors in terms of their 
angular momentum and how the orbital angular momentum 
is re-distributed within the remnant. Ideally one would fol- 
low individual galaxies and calculate the change in \n due 
to mergers and star formation in discs, the latter generally 
contributing to an increase in \r. 

To assess the impact of mergers on the evolution of \r 
one can to first order use the analogy to da rk matter haloes 
and their spin parameter A (jEmsellem et al.ii2007i '). Follow- 
ing the merging history of dark matter haloes, and under 
the assumption of random orbits and complete transfer of 
orbital angular momentum to the remnant, one recovers a 
log-normal distribution f or A consistent wi th results from 
cosmological simulations (jMaller et al.ll2002l ). The distribu- 
tion of A i n dark matter ha loes shows only a weak mass de- 
pendence IjBett et al.ll2007l ). in stark contrast to the strong 
mass dependence of \r in ETGs implying that whatever 
drives the trend in \r with mass is not just simple trans- 
fer of orbital angular momentum. While individual mergers 
might burst or lower \r occasionally based on the orbit 
(Paper VI), continuous merging on random orbits will drive 
\r to a mass i ndependent d i stribu tion at a median value 
of \r ~ 0.038 (|Maller et all 120021 '). Here we used the re- 
lation between the sp in parameter A and Ah presented in 
IjEmsellem et aL 2007) and the median spin parameter of 
dark haloes in lBett et al.l ((20071). Such low \r values would 
classify the majority of ETGs as slow rotators, thus addi- 
tional physical processes need to be important. 

Detailed kinematic analyses of ATLAS'^^galaxies by Pa- 
per II shows, that a large fraction of fast rotators have disc- 
like regular rotation patterns and bars, suggesting that the 
stellar body is hosting a disc component. Stellar discs in 
terms of \r would all be classified as a fast rotating compo- 
nent (Paper VI). If mergers indeed on average lead to low 
Aij, the presence of a disc can be all that makes the dif- 
ference between fast and slow rotators. The process of disc 
building (increasing \r) is in direct competition with the 
process of disc destruction by mergers (lowering \r). 

Motivated by above arguments we introduce a simple 
model based on the mass fraction in discs to investigate the 
origin of fast and slow rotating ETGs. In what follows we 
choose the fraction of fast-rotators as a function of lumi- 
nosity as our key-observable, and gauge our models against 
it. 



4.1 Fast Rotators 

Fast rotators are the majority of early-type galaxies and 
dominate the population at Mk ;i — 24 ( see Fig. [5]). De- 
tailed analyses of several fast rotators by iKrainovic et al.l 
(200q, [201l|) reveal the existence of disc-like features in 
earl y- type host galaxi es in agreement with earlier studies 
of lRix fc White! (|l99CI ). We here assume in agreement with 
Paper VII that such stellar discs are present in all fast ro- 
tators and use their existence as our criterion to distinguish 
between fast and slow rotators. We introduce a free model 
parameter ffast, that gives the lower limit to the disc mass 
in fast rotators. Every early-type model galaxy with a disc 
fraction Mdiac/^'h < ffast is classified as a slow rotator. 
We calculate ii'— band magnitudes for our model galaxies 
using their star formati on history and the stellar popula- 
tion synthesis models of lBruzual fc Charloj (|2003l ). In Fig. 
[5]we show the fraction of fast rotators in the overall popula- 
tion of early-type galax;ies as a function of their rest-frame 
ii'-band magnitude for the ATLAS'^°sample. The sample 
shows a clear trend toward slow rotators at the luminous 
end and is dominated over the remaining observed range by 
fast rotators (Paper HI). We show the fraction of fast rota- 
tors from our SAM for different values of the free parameter 
ffast- In general we find that values of ffast ~ 0.1 — 0.15 
show good agreement with the observations. We here choose 
ffast =0.1 as our fiducial value, based on matching the ob- 
servational data in Fig. [5] and the additional constraint of 
being close the number density of fast and slow rotators from 
the ATLAS^°sample (see Fig.[Tl|- The exact value of ffast 
does not infiuence the shape of the curve shown in Fig. [S] 
it only changes the magnitude at which fast rotators start 
to dominate. The modelled galaxy population of fast and 
slow rotators is in good agreements with the observations, 
supporting the idea that stellar discs are the dividing factor 
between fast and slow rotators. 

The formation path toward fast rotating ETGs in our 
model is either by rebuilding a stellar disc around a spheroid, 
the growth of a spheroid via the destruction of existing discs 
in late-type galaxies or the exhaustion of gas in late-type 
galaxies. Slow rotators on the other hand have two possi- 
ble formation paths: major mergers or repeated dry minor 
mergers. In what follows we will investigate the importance 
of these different paths by changing the model prescription 
within our fiducial SAM. We will focus in particular on the 
importance of gas cooling, the destruction of stellar discs 
in major mergers, star burst during major mergers and the 
effects of minor mergers. 



4. 1.1 Impact of Gas Cooling 

We first change our fiducial model in a way that highlights 
the importance of gas cooling for the formation of fast rota- 
tors. We will call it max-merg in the following (label MM in 
Fig. [6]& Fig. (Tjl. In this model every binary major merger 
with R ^ 3.5 results in a slow rotator, and all the progenitor 
disc mass is added to the bulge, as well as all progenitor gas 
is converted into bulge stars during a star burst. The modi- 
fication to the fiducial model (Eq. [S} is summarised below. 
For major mergers we assume: 



Mbui. 



Mbul,l + Mbul,2 + Mdis.l + Mdis,2 
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+ Mgas,l + Mgas,2 
-^ ''J- gas, rem — U 

and for minor mergers with 10 ^ R > 3.5: 

Mbul,rem = MbuLl + Af*,2 
-'■^-^dis.rem ■'-'^dis.l 



(8) 



(9) 



M„ 



{Mgas,l + Mgas 



In this model the formation of fast rotators via discs is max- 
imally hindered and only allowed via cooling of gas from 
the hot halo and cold gas from gas-rich satellites in minor 
mergers. The solid blue line in Fig. [6] shows the fraction 
of fast rotators that formed via disc growth around bulges. 
We find the expected trend of this mode being more impor- 
tant at lower mass scales, when cooling times are shorter, 
most visible at Mk ~ — 24. However, shorter cooling times 
also mean that more gas will have cooled down before a 
merger takes place, and thus less material is available to 
re- grow a disc. For massive galaxies one can use a rough es- 
timate to evaluate the role of cooling with respect to merg- 
ers. Observations of the merger rate in the nearby universe 
find on average O.OOOlGyr"^ Mpc ~'^ mergers for m assive 
galaxies with M, > 2.5 x 10^° Mg (|jogee et al.ll2009l 1. For 
the ATLAS^^early-type sample this suggests less than one 
merger within the last Gyr per gala xy. Cooling time scales 
on th e other hand are around 1 Gyr (jNeistein fc WeinmannI 
[2OI0I), suggesting that cooling processes and associated star 
formation in disc-like components are not negligible in be- 
tween merger events. However, cooling by itself is not suffi- 
cient to produce enough fast rotators by means of transform- 
ing bulge dominated major merger remnants as evidenced 
by the large drop in the fast rotator fraction at Mk ~ — 24. 
The more important contribution to the fraction of fast ro- 
tators comes from gas-poor low B/T satellite galaxies liv- 
ing in dense environments (see Fig. [T]). On average, fast 
rotators that are 'dried-up' late-type galaxies in the range 
—22 > Mk > —25 are ~ 3 times more frequent than 're- 
grown disc' fast rotators showing that disc re-growth after a 
major merger is a less important channel for the formation 
of fast rotators than possible environmental effects leading 
to the exhaustion of gas in late-type galaxies. 



4.1.2 Gas-rich Mergers 

Within our fiducial merger model stellar discs have two dis- 
tinct origins. They are either the remains of a stellar disc 
that survived a merger or are re-grown from gas that set- 
tled into a disc after a merger. The gas for the latter process 
can be part of the progenitor discs that survived the merger 
or gas that is radiatively cooling and settling onto a disc. 
The importance of the latter process has been addressed in 
the last section and we will turn now to the survival of gas 
in progenitor discs. Numerical simulations show that dur- 
ing mergers of galaxies with high gas fractions, progenitor 
discs are ab le to survive and to s tart quickly rebuildin g stel- 
lar discs CR obertson et al.ll2006l : lHopkins eit al.ll2009al '). Our 
fiducial merger model in Eq. [5] does take the results from 
such simulations into account via the function f burst. In the 
max-merg model above we assumed f burst = 1. To estimate 
the importance that the survival of gaseous disc during ma- 
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Figure 5. The fraction of fast rotators in the overall population of 
early- type galaxies in the ATLAS''^ sample (filled symbols) (Pa- 
per III). The SAM outputs for different values of ffast are shown 
as different coloured solid lines. Throughout this paper we will 
use as our fiducial value for the minimum disc fraction in fast 
rotators ffast =0.1. 



jor mergers has on rebuilding stellar discs we modify the 
max-merg model by using fturst fronr Eq. [5] The survival 
of gaseous disc material during mergers is mainly affect- 
ing the fraction of fast rotators between Mk ~ —23 and 
Mk ~ —24, with only a mild increase in the total number 
of fast rotators by a factor of 1.6 (see Fig. [6]& Fig[7]dotted 
line, label MM_G). From this we conclude, that in our model 
the formation of fast rotators via the growth of stellar discs 
from cold disc gas that is surviving major mergers is less 
important than via radiative cooling of gas. 

4.1.3 Survival of Stellar Discs 

Detailed numerical simulations of binary galaxy mergers 
with and without gas show, that the mass ratio of the 
merging partners plays a crucial role in erasin g the mem- 
ory of the progenitor's stellar disc (e.g. J 3arncs fc Hernguistl 
ll992l : lNaab fc BurkertI [20031 : iHopkinset al.ll2009al ). Within 
our fiducial model we always destroy a fraction of the stel- 
lar disc that is given by min[M£ji,i, M,,2]. Effectively, our 
model states that ten 10 : 1 merge rs have the same effect 
on a stellar disc as one 1 : 1 merger (jBournaud e t al. 2007j). 
To see how the survival of a stellar disc will influence the 
fraction of fast rotators we modify the max-merg model to 
include the dependence on the mass ratio as laid out in 
our fiducial merger model Eq. O Please note that we as- 
sunre that only in major mergers, R < 3.5, the disc of the 
host galaxy will be affected. Minor mergers are neglected 
in terms of their impact onto the stellar disc, and we will 
come back to it in the next section. As can be seen from 
the dashed line labelled MM.D.MAJ in Fig.[6l the fraction 
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Figure 6. The fraction of fast rotators in the overall population of 
early- type galaxies in the ATLAS^'' sample (filled symbols) (Pa- 
per III). The SAM outputs for different model assumptions are 
shown by the various lines in the figure (see text for details on the 
models). For all models we use ffast = 0.1. The blue, green and 
yellow lines show models including the effects of major mergers 
only, while the red line shows the outcome for a model including 
minor mergers. Results in this figure are shown for fast rotators 
that do not have B/T < 0.5 (- low B/T). The survival of stellar 
and gaseous discs during major mergers (models MM_D_MAJ and 
MM_G, respectively) plays mostly a role for fast rotators around 
Mn ~ —23, while the overall population of fast rotators is domi- 
nated by regrowth of a disc from gas cooling and disc destruction 
in minor mergers. 



of fast rotators is mainly increased for luminous galaxies 
around Mk ~ — 23. This increase is modest with respect to 
the max-merg model and of the same magnitude as in the 
case of including stellar disc survival during mergers. Un- 
equal mass major mergers dominate th e overall number of 
mergers (e.g. lKhochfar fc Burkertll200ll ). thus the survival of 
stellar discs is elevated with respect to the max-merg model. 
Again we find that disc re-growth via cooling is more impor- 
tant for the formation of fast rotators after a major merger 
than the survival of pre-existing properties during this major 
merger, in this case the stellar disc of the main progenitor. 
Our results here and in the previous section show that 
surviving discs, stellar and gaseous, are important in allow- 
ing to maintain and re-grow substantial discs from cooling 
of gas in ETGs mostly around Mk ~ —23, and that the ma- 
jority of fast rotators is likely to originate from a different 
formation path within our model which. 



4-1-4 Minor Mergers 

The frequency of minor mergers an ETG encounters during 
its e volution by far exceeds the number of major mergers 
(e.g. lKhochfar fc Silkll2006al : lKavirai et al.ll2009l : lNaab etall 



Figure 7. Same as Fig. [6] but now including the contribution 
from gas-poor low B/T galaxies {+ low B/T). Such galaxies play 
an important role in major merger only models. However, in the 
model including the effect of minor mergers (MM_D_MIN) the 
fraction of fast rotators originated from drying out of late-type 
galaxies in high density environments is significantly lowered. 



12003 : IShankar et al.ll2010l ). The importance of minor merg- 
ers for the formation of bulges in late-type g a laxies has been 
recent ly pointed out (jWeinzirl et al.l l2009l : I Hopkins et al.l 
l2009bl ) and we here investigate their importance in terms 
of fast and slow rotating ETGs. As discussed above, minor 
mergers are able to disrupt parts of th e stellar disc of the ir 
host and to trigger a small star burst (IPeirani et al.ll201Gl ). 
We here modify the max-merg model to include the disrup- 
tion of stellar discs in minor mergers 10 > 7? > 3.5, but 
do not allow for any star burst during it. The solid red line 
labelled MM.D.MIN in Fig. \6\ shows the outcome of this 
model. It is clearly visible how the fraction of fast rotators 
is significantly increased at Mk ~ — 25, with respect to the 
max-merg model. Our results suggest that minor mergers 
are efficient means in gradually transforming disc dominated 
systems into systems with 0.5 < B/T < 0.9, corresponding 
to the B/T range of fast rotating ETGs. As we will show 
below, the continued impact of such minor mergers is lead- 
ing to even higher B/T and the formation of slow rotating 
ETGs. By comparing Fig. [6] and [7] we find that the fraction 
of fast rotators with such an origin is much larger than the 
one for fast rotators that are dried up late-type galaxies. 

4.2 Slow rotators 

Mergers play by construction an important role in the for- 
mation of ETGs in our model. Our assumption of equal mass 
mergers predominantly resulting in slow rotators is clearly 
an oversimplification of what is happening in reality (Paper 
VI). In general, the formation of a slow rotator requires that 
the orbital angular momentum and the spin of the galaxies 
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compensate each other such that only httle angular momen- 
tum is left in the remnant, a situation that requires fine- 
tuning, given that the ratio of orbital angular mome ntum 
to spin is on average larg er than one ( Khochfar fc BurkertI 
l2006d : iMaller et al.1 120021) and that slow rotators after one 
additional major merger tend to spin up to become fast 
rotators again (Paper VI). The question thus arises how 
slow rotators might form in the first place, if even under 
favourable conditions, such as equa l mass mergers, t heir for- 
mation/survival chances are small IjBois et al.ll20ld l. As we 
pointed out earlier, minor mergers are very frequent and 
can contribute a significant fraction of th e remnant mass. In 
an earlier study iKhochfar fc Silkl (|2006al ) showed that up to 
80% of the bulge mass can come from outside the galaxy via 
m ergers, estimates in agreement with numerical simulations 
of lOser et al.l (|201G| ) . If this is indeed the case then slow 
rotators are the end-stage of constant impact by satellite 
mergers, falling in from random directions and helping to 
reduce the angular momentum of the host on average. Since 
the cooling of gas has mostly ceased in high mass galaxies 
living in haloes > MoM.crit, disc formation will be termi- 
nated as well, and thus the conversion into fast rotators. To 
highlight the relative importance of minor and major merg- 
ers we show in Fig. [S]the fraction of minor mergers that the 
most massive progenitors of a present-day slow rotators ex- 
perienced during their evolution. We here define major and 
minor mergers as i? ^ 4 and 20 ^ 7? > 4, respectively, to 
increase the statistical significance. There is a large scatter 
in the fraction of minor mergers that slow rotators experi- 
enced during their evolution and no significant correlation 
with galaxy mass. The median fraction of minor mergers 
that the most massive progenitor of a slow rotators experi- 
ences is ~ 0.7 corresponding to ~ 2.3 more minor than ma- 
jor mergers. The implication of these results is that massive 
slow rotators had most likely a major merger event in their 
past, but got constantly hit by minor mergers, that erased 
any progenitor disc memory and additionally hindered the 
rebuilding of a sufficiently large stellar disc. Looking at the 
individual cases of slow rotators in our model, we find that 
on average slow rotators with high fractions of minor merg- 
ers went through an episode in their past during which they 
were classified as fast rotators, while those slow rotators with 
low minor merger fractions were always close to the slow ro- 
tator regime. 



5 FAST AND SLOW ROTATOR 
DEMOGRAPHICS 

The big advantage of the ATLAS^^sample is its complete- 
ness, which allows us to directly compare our model pop- 
ulation with the observed one. As we have shown in Fig. 
[5] our fiducial model is able to reproduce the relative frac- 
tion of fast and slow rotators as a function of magnitude. 
We find in addition, that the number density of ETGs with 
M* > 6 X 10^ M0 is 1.6 X 10"^ Mpc"^, comparing well 
to the ATLAS^one of 2.2 x 10"^ Mpc'^ (Paper I). Ap- 
proximately 20% of these ETGs are slow rotators, in fair 
agreement with the 14 ± 2% found in the ATLAS^^sample 
(Paper III). Having established that our model reproduces 
the basic statistics of the ATLAS^°sample we can now ad- 
dress the question of what the deciding factors are in the 
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Figure 8. The fraction of minor mergers N minor / {J^ minor + 
Nmajor) that the most massive progenitor branch of a present- 
day slow rotator experienced during its evolution. For the sake of 
better statistics we here define major merger as R ^ 4 and minor 
mergers as 20 ^ i? > 4. The histogram shows the distribution 
of the minor merger fraction which has a median value of 0.73 
(dashed Une) . Minor mergers clearly dominate the history of slow 
rotators. However, we do not find a strong correlation with stellar 
mass. Slow rotators with very few or no major mergers have more 
minor mergers than slow rotators with some major mergers. The 
evolution paths of these two classes of slow rotators is very differ- 
ent with the former going through a fast rotator phase, while the 
latter always stays close to the slow rotator regime. The vertical 
dotted line shows the completeness limit of the ATLAS'^^sample 
(Paper I) 



history of an ETG that make it become a fast or slow ro- 
tator within our model. As we argued above we take the 
disc fraction in ETGs to first order as a proxy of Ah, thus 
the balance between disc destruction and formation plays a 
crucial role. 

We start by looking at the fraction of accreted stars 
that end up in ETGs as a function of their stellar mass. 
By definition these are all stars that a satellite galaxy con- 
tributes to the primary galaxy during a merger and directly 
proportional to the amount of stellar disc that might be 
destroyed in the primary (see Eq. O, suggesting that fast 
rotators should have on average Macc/M* ^ 0.5 depending 
on the amount of gas and stars in the progenitor discs. As 
demonstrated in Fig. [51 the majority of fast rotators has 
indeed Macc/M* ^ 0.5 and the amount of accreted ma- 
terial in slow rotators can reach up to 90% of their total 
mass . These est imates are in agreement with earlier results 
of Khochfar fc S ilk (2006a). When dividing the ETG sample 
into fast and slow rotators a clear separation appears sup- 
porting the observational classification into slow and fast 
rotators. Slow rotators have on average larger fractions of 
accreted material than fast rotators and show a very steep 
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Figure 9. The fraction of accreted stars from satellite galaxies 
onto ETGs as a function of the present-day galaxy mass. Massive 
galaxies accrete the majority of their stellar material. This can be 
up to 90% of the total mass in the most massive galaxies. The 
population of fast and slow rotators show a clear distinction in this 
figure, with the latter (red filled circles) showing Macc/Mt ^ 0.5. 
Slow rotators show a strong trend with mass due to the cooling 
shut-off in haloes more massive than MoM,crit- Fast rotators with 
low B/T ratios are dominantly satellite galaxies in high density 
environments that got stripped of their hot haloes and ceased to 
form stars (light blue filled circles). Please note that we do not 
show satellite galaxies below the ATLAS^'^ survey mass limit of 
~ 6 X 10^ Mq (dashed line), but show central galaxies all the way 
down to 10^ M(?) to make the trend with mass more visible. 



relation between the accreted mass fraction and their total 
stellar mass. 

Fast rotators in our model fall into two different classes, 
those with 0.5 < B/T < ffast (~ 78%) and those with 
B/T ^ 0.5 and low gas fractions (~ 22%). The high B/T 
population shows a clear separation from the other fast rota- 
tors, having higher fractions of accreted material at a given 
stellar mass. Also clearly visible is a moderate correlation 
with stellar mass toward higher accretion fractions. The low 
B/T population on the other hand shows a large scatter 
and stays generally below accretion fractions of 0.3. These 
two populations in terms of their morphology will not nec- 
essarily show different properties other than their B/T ra- 
tio, and thus would to first order look like members of one 
homogenous class of fast rotators. However, detailed anal- 
yses of their stellar population, in terms of age, metallicity 
and detailed structure might reveal differences, due to their 
different formation paths and/or the environment in which 
they evolved. As we showed earlier, the most dominant for- 
mation channel for fast rotators is via the destruction of 
existing stellar discs in minor mergers resulting in high B/T 
fast rotators and possibly thickened stellar discs in contrast 
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Figure 10. The average number of minor (solid) and major 
(dashed) mergers that the most massive progenitors of present 
day ETGs during their evolution experience. The number of merg- 
ers rises with the present-day stellar mass of ETGs. Both pop- 
ulations of ETGs show similar numbers of minor mergers as a 
function of stellar mass. However, slow rotators have on average 
more major mergers than fast rotators. The model estimates on 
the number of mergers is in ve r y goo d agreement with recent 
estimates of lConselice fc Arnold! l|2009l l 



to the low B/T population of dried-up late-type galaxies 
with thin stellar discs. 

The trend seen in Fig. [9] can be understood by consid- 
ering the growrth of galaxies. The ratio Macc/Mt increases 
during dissipationless mergers and is lowered during dissi- 
pational events, like cooling and associated star formation. 
The modest stellar mass dependence seen for high B/T fast 
rotators is mainly driven by the slowing down of cooling 
in massive haloes. Once host haloes reach masses above 
MDM,crit cooling Completely stops and galaxies grow mostly 
via dissipationless mergers (JKhochfar fc Silkl |2009a|) caus- 
ing the strong dependence of Macc/M, on stellar mass seen 
in slow rotators. The low B/T fast rotators resemble the 
flattened fast rotators in the ATLAS sample and are the 
result of cold star forming gas running out in either satel- 
lites galaxies in dense environments (~ 17% of fast rota- 
tors) (light blue filled circles in Fig. O or central galaxies in 
haloes > MDM,crit that did not experience many mergers 
during their evolution (~ 5%) (green filled circles). Another 
20% of fast rotators have high B/T ratios and are satellites 
in dense regions giving a total of ~ 37% of fast rotators be- 
ing satellite galaxies subject to environmental effects (Paper 
I). The majority of ~ 60% however, are central galaxies. In 
contrast, slow rotators are in over 85% of the cases central 
galajcies. 

Besides differences in the fraction of accreted material 
it is important to ask whether the accreted material is pro- 
vided in the same way for fast and slow rotators. In Fig. 
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[To] we show the average number of major and minor merg- 
ers the most massive progenitor of present-day fast and slow 
rotators experienced. Not surprisingly, we find that the num- 
ber of mergers is an increasing function of stellar mass and 
that the average number of minor mergers is similar for fast 
and slow rotators of the same mass. The number of major 
mergers however, is larger for slow rotators. The average 
number of major mergers is between 1 and 2, which implies 
that they encountered at least one major interaction during 
their formation which would leave its imprint. 

In summary, fast and slow rotators mark the transition 
point in the ability for gas to cool and form stars in discs. 
Slow rotators have accreted more stars via minor and major 
mergers, reducing their disc fraction and hence Xr. In addi- 
tion they have more major mergers, at least one on average 
during their evolution that could leave an imprint in their 
stellar body. 
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6 THE REDSHIFT EVOLUTION OF FAST 
AND SLOW ROTATORS 

As we laid out earlier, in our model the ability to cool gas 
and the fraction of accreted stars via mergers are important 
to predict the fraction of fast and slow rotators. The impact 
of both these physical processes is changing going to higher 
reds hifts. The merger rate of galaxies peaks around z ~ 2— 3 
fe.g. IConselice et aLll2003l ), and the gas accretion rate als o 
does increase towards higher redshifts (jDekel et al.l 12003 ). 
In Figure [11] we present the number density of fast and slow 
rotating ETGs as a function of redshift for three different 
ranges in stellar mass (9.8 < logM* < 10; 10 > logM, < 11 
and logM* > 11). The horizontal line shows the number 
density of early-type galaxies in the ATLAS'^'^ volume. Going 
to higher redshift the number density drops as expected due 
to the hierarchical build-up of massive galaxies that takes 
place, and the fact that a larger fraction of massive galaxies 
are classified as gas-rich late-type galaxies. Independent of 
the overall decline, the number density of slow rotators in 
all mass bins in general declines stronger with respect to 
the one of fast rotators. Many of the slow rotators at high-z 
are remnants of very recent equal mass mergers, in contrast 
to the low-z slow rotators, that have had major mergers 
in their past but are dominated by minor mergers during 
the later stages of their evolution. As shown by Paper VI 
equal mass mergers that take place under special merger 
orbits can result in fast rotators. The number densities for 
slow rotators at high-z that we present are therefore upper 
limits. Independently, our results predict that beyond z ~ 2 
it will be hard to find slow rotators and that the population 
of ETGs at 2 ^ 2 should be dominated by fast rotators, 
which are more than one magnitude more frequent than slow 
rotators at M« > 10^" M0. Taking the observational results 
on the evolutio n of the number-den sity of massive ETGs at 
face value (e.g. JFerreras et al.ll20C)9r ) our results suggest that 
basically all massive slow rotators were fast rotators at some 
point during their evolution. 



Figure 12. The FOF-mass function of dark matter haloes 
in a cosmological N-body simulation (red line) of volume (100 
Mpc/h)'', and in 50 spherical sub- volumes of the size of the 
ATLAS'^'^volume, placed randomly within the large simulation. 
The effect of cosmic variance is clearly visible at the massive end 
of the mass function. 



7 COSMIC VARIANCE 

The ATLAS^^— survey only spans a limited volume of the 
local universe within a radius of 42 Mpc. It is thus valid 
to ask whether cosmic variance has a significant impact on 
the results of the survey. We try to address this question 
by using our fiducial model and investigating the effects of 
cosmic variance using a large scale N-body simulation. 



7.0.1 Variations in the Dark Matter Mass Function 

Depending on the location of a fixed volume within the 
Universe, the mass function of dark matter host haloes 
and associated galaxies will change. We estimate this ef- 
fect, by using a cosmological N-body simulation of a (100 
Mpc//i)^ volume. The particle resolution of the simulation 
is 10^ M0//1, and we generate the z — Q mass function 
for haloes using the friends-of-fiends (FOF) algorithm. Fig- 
ure [12] shows the mass function of the whole box (red line) 
and that of 50 randomly placed spherical sub-volumes of 
the size of the ATLAS'^° volume. The mass functions agree 
well at low masses, but starts deviating at large masses 
due to cosmic variance. To predict the fraction of fast and 
slow rotators within each of these sub-volumes we take the 
FOF-mass function and ap ply the merger-tree algorithm of 
ISomerville fc KolattI ([1999|) in combination with our fidu- 
cial model. We do not construct the merger trees from the 
simulation, because we do want to also investigate what the 
impact of randomly generated merger histories is, choosing 
haloes of the same mass (see § 17.0. 2p . 

The fraction of fast rotators in Fig. [13] shows fiuctua- 
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Figure 11. Number density of fast and slow rotators as a function of redshift. We apply three different stellar mass bins to highlight the 
mass dependence. The long dashed lines show the number-density of slow (red), fast (blue) and all (black) ATLAS^'^ galaxies at 2 = 0, 
respectively. Our model predicts a strong decline of slow rotators even for the highest mass bin and predicts that fast rotators should be 
more than an order of magnitude more frequent at ^ > 2 for galaxies with M« > 10^*^ Mq. The high-z number densities of slow rotators 
are upper limits since they are the results of recent equal mass mergers that under certain merger orbits result in fast rotators (Paper 
VI). 



tions at Mk ~ — 25, which are well within the error-bars of 
the observed sample. Only for more luminous galaxies the 
fraction of fast rotators is heavily affected by cosmic vari- 
ance. However, the fraction of fast rotators in the highest 
luminosity bin is on average still below 30%. In terms of the 
ATLAS''°sample our results suggest that the general trend 
of a declining fraction of fast rotators as a function of lumi- 
nosity is robust and not affected by cosmic variance. 



ble to the one seen for different mass functions and suggests 
that part of it might be initially due to the varying merger 
histories and only to second order by the different number 
densities of haloes in a given volume. Varying merging his- 
tories has most impact at the extreme high-mass tail, where 
the relatively small number of major mergers can be strongly 
influenced by a few more major mergers during the history 
of a galaxy. 



7.0.2 Variations in Merger Histories 

Besides changes to the dark matter mass function in a given 
volume the individual merging history of the dark matter 
haloes might be different as well, and have a significant im- 
pact on the population of early-type galaxies. We investigate 
the impact of the merger histories by generating 30 random 
merging histories for the dark matter mass function drawn 
from one of the sub- volumes. Again we apply our fiducial 
model to predict the fraction of fast and slow rotators. The 
scatter about the median value is well within the error bars 
of the observations (Fig.ll4p. and large deviations only occur 
at Mk /$ — 25. The magnitude of the scatter is compara- 



8 DISCUSSION & SUMMARY 

In this study we have used a semi-analytic modelling ap- 
proach that by construction reproduces first order key 
observations, such as the mass function, the trend in 
bulge-to-total stellar mass ratio as a function of stel- 
lar mass and the gas fraction in late-type galaxies. On 
top of this model we have made a selection of ETGs 
that tries to be as close to the ATLAS^ ° (Paper I) selec- 
tion as possible. While historically B /T (|Kauffmann et al.l 
1 19991 : ISimien fc de Vaucouleurdl 19861 ) has been the preferred 
choice, in this study we extended this criterion to take into 
account the revised definition of ETGs proposed by the 
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Figure 13. Same as Fig. [5] for 50 random spherical sub- 
volumes within a large scale cosmological N-body simulation. 
The cosmological simulation itself has a volume of (100 Mpc/h)'' 
and each sub-volume has a volume comparable to that covered 
by the ATLAS^° sample, i.e. Va3D = 10^ Mpc^. Our results 
show that cosmic variance should not play a large role for the 
ATLAS'^° sample at Mk > —25. Only the fraction of fast rotators 
at the very luminous end might be affected by cosmic variance. 
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Figure 14. Same as Fig. 1131 but for 30 random merging trees 
using the same mass function of haloes at 2 = 0. Below Mx ~ 
—25 results are robust and within the errors. The fraction of fast 
rotators at very low and very high luminosity is heavily affected 
by the individual merging histories of dark matter haloes. 



ATLAS"^°survey (Paper I). We introduce the additional con- 
dition of low cold gas fractions in disc dominated galaxies 
to account for the population of flattened gas-poor fast ro- 
tating ETGs in the ATLAS^°sample. Support for this ap- 
proach comes from the observed distribution of gas fractions 
in late-t ype and early-typ e galaxies in the ATLAS^'^parent 
sample IjSerra et al.l [201l]) , and from the fact that binary 
mergers in general do not result in fast rotators as flat as 
observed (Paper VI). Anyway, the fraction of model fast 
rotators with low B/T and low gas fractions is only 22% 
and is mainly affecting the low mass end, where fast rota- 
tors dominate in any case. We thus expect this ETG se- 
lection not to change any trend seen between the fraction 
of fast and slow rotators at high masses. Based on these 
ETG selection criteria our model predicts a number den- 
sity of ETGs in good agreement with the ATLAS^'^parent 
sample. We did not apply any 'fine tuning' to match the 
number densities of ETGs, and these are just the outcome 
of the selection procedure that we applied. It is interesting 
to note, that more than 80% of ETGs in our sample have 
B/T values consistent with the classical selection of E and 
SO galaxies based on the Hubble T-type an d bulge-to-total 
light ratios (|Simien fc de VaucouleursHigSq ). Similar agree- 
ment has been found between the T-type based morpho- 
logical definition and the ATLAS^°ETG selection (Paper I) 
further supporting that our model ETG selection is close to 
the ATLAS^°one. 

In our model we grow bulges only via mergers and we 
neglect internal processes like disc instabilities that have 
been suggested as additional ways to grow bulges (e.g. 
IParrv et al.] l2009l : iDekel et all |2009| ). These processes are 
most efficient in gas-rich massive discs and would result in an 
increased fraction of fast rotators with respect to our fiducial 
model with a range in B/T ratios depending on when the 
disc becomes stable. In such models clumps form in the disc 
and migrate due to dynamical friction into a central bulge. 
It is however, not clear at this point whether such clumps 
survive long enough to rea ch the centre or ge t disrupted by 
feedback from supernovae IjGenel et al.ll201Cl ). 

Another important ingredient in our model is the abil- 
ity to cool gas and re-grow discs. While SAMs in general 
reproduce the average cooling behaviour of gas in SPH sim- 
ulations there ar e difference when comparing on an object- 
to-object basis (|Saro et al.l [2010). Even between individ- 
ual SAMs deviations in the calculated cooling rates exist 
l|De Lucia et al.ll201(J ). In general, the fraction of massive 
slow rotators is not affected by changes in the cooling rate, 
due to the constraint that the galaxy mass function needs to 
be reproduced. This is usually achieved by feedback imple- 
mentations in the SAMs, which regulate the cooling of gas 
and the overall star formation. SPH simulations show that 
feedback plays indeed a key role in regulating gas accre- 
tion rates and star formation (e.g. Ivan de Voort et al.ll2010l : 
ISchave et al.ll2010l ). In our model we 'tune' our feedback 
efficiency and hence our cooling rate in such a way that 
we reproduce by construction the following observables: the 
stellar mass function, the gas fraction in late-type galaxies 
and the cooling rate in galaxies at high z. Using this ap- 
proach we try to minimise the freedom in our cooling rates. 

Besides cooling of gas, mass loss from already exist- 
ing stars can provide significant material for the growth of 
stellar discs JMartig fc BournaudI [20101 : lAgertz et al.ll2010l : 
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iKimm et al.ll2011h . We do not take this into account in our 
model, but note here that s uch mass loss is likely to dec rease 
the B/T ratio of galaxies. i Martig fc BournaudI ()2010l ) find 
that the Hubble T-type of a galaxy can change via this pro- 
cess by 2 to 3 units, affecting the fraction of fast rotators. 

The merger rate of galaxies is a key aspect in our mod- 
elling approach. Being able to reproduce the stellar mass 
function, we populate dark matter haloes on average with 
galaxies of the right mass. Thus when dark matter haloes 
merge we expect to merge the 'right' galaxies as well. In 
fact, the derived merger rates and average number of merg- 
ers from our model are in very good agreement with obser- 
vations (jJogee et al.ll2009l : IConselice fc Arnold! l2009l ). 

Within the ETGs selection we distinguish between fast 
and slow rotators based on the fraction of stellar disc still 
present in the main body. The main reasoning behind this 
approach is the observational fact that many fast rotators 
show regular rotation patterns and signs of bars (Paper II), 
which indicate the presence of a disc-like component (Paper 
VII). In addition disc galaxies show values of Xr > 0.5 sim- 
ilar to fast rotating ETGs in the ATLAS^°sample (Paper 
VI). In Paper VI the results of binary mergers have been 
extensively compared to the structure of galaxies in the 
ATLAS"^°sample showing that binary disc mergers repro- 
duce best the properties of fast rotators with Xr ~ 0.25—0.5, 
while fast rotators with Xr < 0.25 are mostly reproduced 
in re-mergers of disc merger remnants, and that the ma- 
jority of slow rotators cannot have formed from a binary 
major merger or re-merger. They in general have a much 
more com plex formation hi story that will be investigated in 
iNaab et a l. (2011) and Boi s et al.l (|2011bl ). Paper VI shows 
that major mergers efficiently reduce Xr in the central re- 
gions. Minor mergers on the contrary are expected to re- 
duce Xr in the outer regions. These results suggest that 
within the fast rotator population a sequence of mergers 
will preferentially lead to a reduction of Xr. The way this 
can be naturally achieved is by the destruction of progeni- 
tor disc components via violent relax;ation during mergers. 
However, discs can re-build by gas accretion and star forma- 
tion increasing A^, and it is therefore important to model 
the complete merging history of a present-day ETC to de- 
cide on its status in terms of fast or slow rotation. We here 
do not attempt to model the detailed structures of fast and 
slow rotating ETGs but try to give answers on the average 
formation path of these galaxies. To achieve this we intro- 
duce a free model parameter, the minimum disc fraction 
required in fast rotators to divide our model ETG sample 
into fast and slow rotators. This parameter has been set to 
ffast =0.1 by matching the fraction of fast and slow rota- 
tors as a function of magnitude and their number density 
at z = 0, and serves the dual purpose of a free model pa- 
rameter as well as a prediction of the expected disc fraction 
in fast rotators that can be compared to observations, once 
proper bulge-disc decompositions of ATLAS'^^galaxies have 
been performed (ATLAS^^in prep.). Effectively our model 
assumes that Xr is to first order a proxy for the disc fraction 
in ETGs. 

The modelled fast and slow rotators show clear distinct 
growth histories, supporting the observationally motivated 
separation into two classes. The two main differences are 
that slow rotators grow by more accretion of stellar mass 



from satellites than fast rotators of the same mass, and that 
they have on average more major mergers. 

The fraction of accreted material depends strongly 
on the competing effects of disc growth in galaxies and 
their destruction during mergers. While fast rotators con- 
tinue accreting small levels of gas and forming stars in 
their recent hi story, consistent with observed star forma- 
tion in ETGs (Schawin ski et al.l 12003 ). slow rotators com- 
pletely stop accreting gas, and only grow via mergers. 
The main mechanism causing cooling to stop in high mass 
galaxies is still debated and could have several different 
origins (see e.g. Croton et al. _ 20061: Catta,neo et al.l 120061 : 



iKhochfar fc Ostrikeill2008l : [Johansson et al.ll2009l ). We here 
adopt an empirical sharp transition in the cooling behaviour 
of gas based on a critical dark halo mass. It is likely that 
such a transition should in reality be smooth based on the 
underlying physics. Introducing a scatter in the critical halo 
mass, to model a smooth transition, is not changing the 
general growth properties of fast and slow rotators though. 
It is therefore safe to assume that whatever physical pro- 
cess causes cooling to stop in massive haloes, as long as it 
scales with galaxy mass, respectively halo mass, it will not 
change the growth history of fast and slow rotators. Com- 
pletely neglectin g any cooling shut-o ff results in too many 
massive galaxies (jCroton et al.ll2006l ). and in particular too 
many massive fast rotators. 

The higher accretion fraction in slow rotators has an- 
other direct implication for the hot halo gas. The po- 
tential energy of infalling satellite galajcies is an effec- 
tive source of gravitationa l heating of the hot halo gas 
(|Khochfar fc Ostriked l2008l : I Johansson et al.ll2009l). Ei ther 
via dynamical friction heating ( El-Zant et al.l [200J) or 
shocks large parts of the potential energy transfers to the 
halo gas. Our model thus naturally predicts to find prefer- 
entially hot X-ray haloes around slow rotat ors, and not fast 
rotators of the same mass (jSarzi et al.ll20ld ). suggesting that 
X-ray haloes are a consequence of the mass assembly of an 
ETG. 

Another important consequence of the higher accretion 
fractions in slow rotators is that even though they have the 
same average number of minor mergers as fast rotators of 
the same mass, the masses of the satellites are higher and 
the relative contribution towards lowering A_h is much more 
important. Higher accretion fraction means higher disc de- 
struction fraction through violent relaxation during the evo- 
lution of a galaxy, and thus lower Xr. We find on average 
less than 10 minor mergers for the most massive progenitors 
of fast and slow rotators (Fig. Ill|) . which allows the mass 
distribution of merging satellites to be different for fast and 
slow rotators. 

The average cold gas fraction of satellite galaxies merg- 
ing with progenitors of present-day slow rotators is lower 
than that of fast rotators. This bias is a direct consequence 
of the fact that the infalling galaxies in slow rotators are 
more massive and thus have lower gas fractions (see Fig. Q, 
which hinders additionally the rebuilding of discs. 

The second main difference in the growth history of fast 
and slow rotators is the number of major mergers. Slow rota- 
tors have on average more than twice as many major mergers 
than fast rotators. However, even the most massive ones do 
not have more than 3 major mergers on average during their 
evolution. We find that most of these major mergers happen 
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at early times between gas-rich progenitors compar able to 
those simulated in Paper VI and IWuvts et al.1 (|2010r ). Such 
gas rich major mergers produce kinematically distinct cores 
(KDC) similar to ones observed in present-day slow rotators 
(Paper VI). However, the remnants from major merger sim- 
ulations tend to be not round enough compared to observa- 
tions. Based on the merging history in our SAM we find that 
many minor mergers follow such a maj or merger event. Mi- 
nor m ergers from random directions (jKhochfar fc BurkertI 
[2006g), if not too compact, get stripped of their stellar ma- 
terial in the outer parts of the host galaxy, providing an 
envelope of stars that can make too flat major merger rem- 
nants with KDCs become more round. Further indirect ev- 
idence for the importance of minor mergers for the growth 
of a host galaxy particularly in its outer l ayers comes from 



the observed size-evolution of ETGs (e.g. Khochfa r fc Silk 



2006b; Truiillo ct al. 2007; Bezanson et al.1 2009: Na ab et al 



20091 : iHopkins et aLlbOlOd ). As shown in Paper VI KDCs are 



very fragile during major mergers. Their observed presence 
in slow rotators prompts the question how they can survive 
till today. Analysing the growth history of slow rotators in 
our model, we find that in most of the cases the last gas-rich 
major merger w as not very recent (z > 1.5) (see however, 
IDuc et al.l ([20111) for the case of a remnant from a possible 
recent gas-rich major merger), and that the number of major 
mergers is just not very high (see Fig. I10|) . both favouring 
the survival of KDCs in slow rotators. In contrast fast ro- 
tators have on average less than one major merger in their 
history, making them unlikely candidates to host large-scale 
KDCs. 

Within our model we predict a strong evolution of the 
slow rotator fraction toward low z. At z ~ 2 the number den- 
sity of massive (A/, > 10^" Mq) fast rotators is more than 
one order of magnitude larger than the slow rotators one. 
This is mainly due to high cooling rates and gas fractions 
in mergers promoting stellar disc building over destruction. 
The evolution toward low redshifts is driven by the inability 
to efficiently re-build stellar disc in massive host galaxies, 
which is another representation of down-sizing in star for- 
mation. 



9 CONCLUSIONS 

In this paper we made use of the completeness of the 
ATLAS''°sample to investigate the origin and formation his- 
tory of fast and slow rotator ETGs within a self-consistent 
cosmological framework using the semi-analytical modelling 
approach. We here present a model in which the difference 
between fast and slow rotator ETGs is purely based on the 
stellar disc fraction found in them, and predict that fast 
rotators have disc fractions Mdisc/Mt, > 0.1. We find that 
slow rotators within an evolving universe mark the tran- 
sition in the ability of galaxies to cool gas and to rebuild 
stellar discs. We find a clear separation in the growth his- 
tory of fast and slow rotator ETGs, supporting the obser- 
vationally motivated distinction. In particular the accreted 
fraction of stars shows a clear distinction between fast and 
slow rotators, with the latter having between 50% — 90% of 
their stellar mass accreted from satellites while the former 
has < 50% accreted. 

Although we find a clear separation into fast and slow 



rotators, we also find that ETGs can switch their state be- 
tween fast and slow rotator and vice versa based on stellar 
disc growth or destruction by mergers, suggesting that fast 
and slow rotator ETGs are transient. These changes how- 
ever, occur predominantly at higher redshifts when cooling 
is more efficient and mergers more frequent. The fraction of 
slow rotators, shows a continued increases with time due to 
the conversion of fast rotators and the inability of gas to 
cool and convert slow rotators back to fast rotators. Mas- 
sive present-day slow rotators therefore, can be viewed as 
the final stage in the evolution of ETGs. 

Future high redshift observations of ETGs will be able 
to reveal any possible evolution in the fraction of fast and 
slow rotators and test the presented model further. 
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